Spectral trends in the Very High Energy blazar sample 

due to EBL absorption 

Bagmeet Behera and Stefan J. Wagner 

Landessternwarte, Universitdt Heidelberg, Konigstuhl, D691 17 Heidelberg, Germany 

Abstract. The absorption of y-rays in the intergalactic medium due to the EBL (Extragalactic Background Light) causes 
the observed blazar spectrum to be fainter and softer than their intrinsic state. It could thus be expected to see an effective 
spectral-softening trend with redshift. No such trend is evident in the sample of VHE blazars currently observed. 

To check which distributions of the properties of the parent blazar population could reproduce the observations, various 
simulations are done. The resulting subsamples that satisfy a generic detection criterion for the current generation of ACTs 
(Atmospheric Cherenkov Telescope) are checked to identify whether any inherent correlations (of spectral properties with 
redshift) are required to explain the current observations. 
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1. INTRODUCTION 

Advances in ground based ACT (Atmospheric 
Cherenkov Telescope) techniques in the last two 
decades has resulted in a dramatic increase in the num- 
ber of extragalactic objects detected in VHE (Very High 
Energy, defined as E> 100 GeV) 7 -rays. The VHE 
blazar sample has increased from just a few, to over 20 
at present. Most of the detected AGN are classified as 
HBLs (High frequency peaked Bl-Lac object); however, 
two LBLs (Low frequency peaked Bl-Lac object), one 
nearby FR-I (Fanaroff Riley type-I) radio galaxy, and 
one FSRQ (Flat Spectrum Radio Quasar) have now been 
detected. The maximum redshift (z) of a source detected 
with high significance is 0.536 for 3C279. However, 
this source has been detected only within a one day 
period, from all observations made public. The redshift 
(confirmed) of the farthest, relatively steady source is 
0.212 for 1ES 1011+496. Thus the AGN zoo is now 
open for study in the highest energy regime. 

VHE emission from all extragalactic sources suffer 
extinction in the intergalactic medium due to the EBL 
(Extragalactic Background Light). This extinction is the 
result of pair-production from photon-photon interac- 
tions (Yvhe + Yebl — * e + + having the most ef- 
fective cross-section for the product of photon ener- 
gies » 1 TeV x 1 eV (i.e. the VHE photon energy taken 
in TeV and the EBL photon energy taken in eV). The 
EBL in the UV-optical and the IR band provide enough 
photon density to get noticeable absorption of source 
photons in energies ^ 100 GeV, and this absorption in- 
creases sharply above a few TeV. The net optical depth 
(t) in this process is naturally dependent on z, giving 
T = t(z,Evhe)- The absorption of y-rays, hence causes 



the observed blazar spectra to be softer than their in- 
trinsic state, with a more pronounced effect for large-z 
sources. 

In this work, the VHE blazar sample is studied for any 
discernable trends in the intrinsic and observed spectral 
features with z. It has been attempted to explain the spec- 
tral trends or the lack thereof, by simulating parent sam- 
ples with various intrinsic properties, that are either de- 
rived from the observed sample, or assumed from certain 
blazar-unification models. 



2. VHE BLAZAR SAMPLE 

The VHE blazar sample as known at the time of this 
symposium is described in Table [TJ Columns 3 and 4 
are the observed photon index and the normalisation at 
1 TeV to a simple power law fit (^ = NQ.fj t E~ T ° bs ), 
respectively. The fit to the observed data was converted 
into an intrinsic spectra by correcting for the attenuation 
due to the EBL; the upper limits given in lfL2ll was 
the template EBL used for this work. The slope of the 
corresponding deabsorbed spectra (calculated by taking 
a straight line fit between 0.2 TeV and 1 TeV) is given 
in column 6, is the intrinsic photon index. The table is 
sorted by redshift. It is clear from this table as well as 
Figure|2]and Figure[3](squares) that there is no evidence 
of a spectral trend with z, for both the intrinsic and the 
observed photon indices. 

This might be the result of the actual distribution of the 
intrinsic properties of the parent blazar sample, which 
this small sample of detected sources correctly reflects. 
However, it might also be due to selection effects, since 
all these sources were carefully selected for pointed ob- 



TABLE 1. The VHE blazar sample, sorted by redshift. The values for the observed photon indices and the 
normalization (to a power law fit) at 1 TeV are from the references in the last column. The Integral flux in column 
5 is calculated from the power law fit with the given r o ^ s and normalization, in columns 3 and 4 respectively. The 
6th and 7th column give the approximate calculated values (corrected for EBL absorption) for the intrinsic T and 
the deabsorbed integral flux between 0.2 TeV to 10 TeV. 
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FIGURE 1. Distribution of the various properties of the VHE blazar sample. Left: The distribution of the intrinsic F (EBL 
absorption corrected) and the Gaussian fit (with mean at 2.2) to it. Right: The EBL corrected flux (0.2 TeV < E < 10 TeV), and the 
powerlaw fit (with slope of 15.481) used for simulations. 



servation with ACT instruments based on certain criteria, 
such as described in [ 171 . and not the result of an unbi- 
ased all sky survey. In this work, a preliminary attempt 
has been made to explore the first of the two possibilities 
mentioned above. 



3. SIMULATING VHE BLAZARS 

Neglecting at first the selection effects involved in choos- 
ing blazars for VHE observations with ACT, various par- 



ent samples are simulated. The properties used to char- 
acterize each parent sample are the distribution of the 
intrinsic photon index, the deabsorbed band flux (as a 
proxy for the intrinsic luminosity), and the redshift distri- 
bution (i.e. the density evolution). The deabsorbed pho- 
ton index (r,„,) distribution is well described by a Gaus- 
sian, Figure [T] left panel, and the band flux (L^y, be- 
tween 0.2 TeV to 10 TeV) is fitted with a power-law in 
N versus Log(L7- e y). Each simulated blazar is uniquely 
represented by a combination of z, r,„ ; and Lr e v, where 
the value of each quantity is obtained from simulations of 
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FIGURE 2. Results of simulation of a 1000 blazars; squares (in both panels) are the true VHE blazars. Left: The intrinsic 
r versus the deabsorbed band luminosity (0.2TeV<E< lOTeV) of the simulated parent population (limited to deabsorbed 
LreV > 10~ 12 /cm 2 / s). Inverted triangles denote the parent sample with the luminosity and T distribution derived from the 16 
VHE sources (Figure [TJ. Circles denote the parent population with the mean intrinsic-r shifted to 2.7. Right: The simulated 
sources that survived the HESS sensitivity cut (368 circles and 309 triangles) on the EBL-absorbed Flux (above 0.2 TeV). There is 
a clear softening trend of the observed T with redshift, in both cases. 



particular functional forms of the various distributions as 
described below. For all sets of simulations an uniform z- 
distribution is taken. Blazars in the simulated samples are 
individually corrected for intergalactic attenuation due to 
the EBL and a detector sensitivity cut (for simplicity only 
a cut to HESS sensitivity is considered here) is applied 
to get the subsample of blazars that could be detected by 
the current generation of ACTs. The resulting subsam- 
ples are then compared to the true VHE blazar sample. 
Three sets of simulations were considered: 

1) For the first set, the functional forms shown in Fig- 
ure [T] were used for the hjev and ^int distributions to 
simulate 1000 artificial blazars (inverted triangles in Fig- 
ure p}, which were considered as the parent sample. The 
absorbed spectra was calculated, using the relevant opti- 
cal depths calculated from the EBL model taken from 
|[T2l . The integral flux in the energy range 0.2 TeV to 
10 TeV of the resulting spectra was found. Sources with 
the integral flux greater 10~ 12 /cm 2 /s were considered to 
be detectable with the present generation of Cherenkov 
telescopes. The fact that no value in column 5 of Table [T] 
is lower than this value, demonstrates that this is indeed 
an accurate description of the detection threshold. Re- 
sults are shown in the right panel of Figure [2] as inverted 
triangles. 

2) The parent blazar sample for the second set was 
chosen with the same Lj e v distribution as above but with 
the mean of r,„, distribution softened to 2.7 instead. The 
other parameters for the Gaussian fit to Yi nt were not 



modified. The same recipe was followed for extracting 
the EBL-absorbed spectra and the get the sub sample of 
sources that can be detected (see Figure[2] circles). 

3) The parent sample for the last set of simulations was 
chosen with the same Lj e v as above but with the mean 
of Tim distribution chosen according to a simple inverse 
linear relation between LagiLjev) and r,„,, representing 
the so called "blazar sequence". According to the blazar 
sequence, 1181 . {19\ and references therein, blazars with 
higher luminosity (over all A, i.e. the bolometric lumi- 
nosity) have the synchrotron peak at lower energies, and 
have a softer spectral index. Higher overall luminosity 
with softer spectral index roughly translates to lower 
band flux in 0.2 TeV to 10 TeV, since brighter sources 
with softer spectra drop sharply at high energies (i.e. 
when extended to VHE) whereas fainter sources with 
harder spectra extend further into VHE before falling off. 
Hence the relation seen in Figure|3] left panel. To get the 
subsample of blazars that can be detected with ACTs the 
same procedure as the above cases was applied. 

4. RESULTS AND DISCUSSION 

The first two simulations show a softening trend of the 
observed photon index with redshift as seen in Figure [2] 
(for both the true distribution as well as a softer range 
of Ti nt ) contrary to the lack of any similar trend in the 
redshift-T relation of the 16 known VHE blazars. The 
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FIGURE 3. The simulation of 1000 blazars (filled circles); squares are the true VHE blazars. Left: The intrinsic F versus band 
luminosity of the parent sample, clearly seen is the inverse relation between Log^^y) and T assumed for this run. Right: The 
simulated sources that survived the HESS sensitivity cut (335 in total) on the EBL absorbed Flux (above 0.2 TeV). There appears 
to be no hint of a spectral softening with z for this simulation. 



third simulation (analogous to the so called "blazar se- 
quence") does not show any such trends. 

It should however be pointed out that the blazar se- 
quence is derived from a highly (observationally) biased 
sample, and might not reflect a true systematic trend in 
the physical properties of AGN. Also, note that the in- 
strumental bias in the sample of the VHE blazars is not 
uniform since the measurements were taken with differ- 
ent instruments having different sensitivities. An addi- 
tional caveat regarding the simulations is that the proxy 
used for luminosity, i.e. the band flux has a z dependence 
and hence is not a completely independent intrinsic pa- 
rameter for the simulated blazars. Simulating the intrin- 
sic band luminosity instead of the band flux produces 
identical results, and will be the parameter used in fu- 
ture simulations to be published elsewhere. Despite these 
shortcomings the results of these preliminary simulations 
are encouraging, as it provides evidence that the presence 
or the lack of a spectral softening trend with z, could be 
the result of a peculiar parent population, though the in- 
strumental bias has to be carefully studied. 

5. CONCLUSION 

We tentatively conclude that the observational bias and 
instrumental-sensitivity limitations of Cherenkov tele- 
scopes, are responsible for cancelling out the expected 
trend of spectral softening (due to EBL attenuation) with 
z, in blazars detected in VHE y-rays. More detailed sim- 
ulations and statistical tests will be done to make quan- 



titative conclusions about the current VHE observations, 
and to make predictions about the trends that might ap- 
pear with a larger VHE blazar sample. 
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